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The paramagnon contribution to the quasiparticle inelastic scattering rate in disordered supercon- 
ductors is presented. Using Anderson's exact eigenstate formalism, it is shown that the scattering 
rate is Stoner enhanced and is further enhanced by the disorder relative to the clean case in a 
manner similar to the disorder enhancement of the long-range Coulomb contribution. The results 
are discussed in connection with the possibility of conventional or unconventional superconductivity 
in the borocarbides. The results are compared to recent tunneling experiments on LuNi2B2C. 
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I. INTRODUCTION 

In the early seventies, certain rare earth ternary com- 
pounds were found to display superconductivity while 
at the same time showed strong tendencies to be mag- 
netic. A large body of theoretical work has been de- 
voted to the interplay of magnetism and superconduc- 
tivity j^]. Recently, there is increasing evidence that 
there is an interplay of magnetism and superconductiv- 
ity in the boro-carbides Q as well as RuSr2GdCu208 
PI . Presently it is widely debated whether these materi- 
als are conventional superconductors with sharply peaked 
density of states (DOS) near the Fermi level (similar to 
some A-15's), or unconventional with ground state pair- 
ing of lower symmetry than the underlying lattice. For 
instance, in i?=Lu, Y i?Ni2B2C, scanning tunneling mi- 
croscopy (STM) has given evidence for conventional BCS 
behavior, albeit with a substantially smeared DOS, with 
a smearing parameter F/A = 0.2 Q. Optical conductiv- 
ity studies also support moderately strong coupled con- 
ventional superconductivity with 2A/rc = 3.9 - 5.2 [|. 
At the same time de Haas- van Alphen |^ , magnetic field 
anisotropy , and electronic Raman scattering exper- 
iments have given evidence for at least very small gaps 
over a portion of the Fermi surface. 

Substitutional or positional disorder has played a cru- 
cial role in determining whether a material is conven- 
tional (i.e., obeys Anderson's theorem) or not, and re- 
cently studies on borocarbides doped with Co have been 
performed Heat capacity and magnetic measure- 

ments on Y(Ni2_xCo2:)B2C have interpreted the 
drop in Tc with increasing Co doping as due to the re- 
duction of the DOS at the Fermi level rather than pair- 
breaking by nonmagnetic impurities. On the other hand, 
Raman measurements on the same systems have shown 
an increase in spectral weight below the gap edge as Co 
is doped in, contrary to conventional BCS behavior ||]. 

However it is well known that conventional supercon- 
ductors which are highly disordered display substantially 
smeared BCS properties which can mimic unconventional 
pairing [ pT| . This can result from vanishing of phase co- 
herence or from the interplay of interactions and disorder. 
The latter is most manifest in the reduction of quasipar- 
ticle (qp) lifetimes. Inelastic collisions broaden qp eigen- 



states and lead to a smearing of activated or threshold 
behavior in single- and two-particle correlation functions, 
measured e.g., by tunneling, optical conductivity, and 
electronic Raman scattering. While the present status of 
the superconducting ground state of the borocarbides re- 
mains unclear, it is of interest to inspect whether strong 
inelastic scattering can modify s-wave properties to the 
point where the ubiquitous exponential behavior of vari- 
ous thermodynamic and transport quantities is obscured. 
For instance, the absence of a coherence peak in NMR 
is usually interpreted as a signal of unconventional elec- 
tronic pairing. However, it is well known that the coher- 
ence peak can be suppressed as a consequence of strong 
inelastic electron-phonon collisions jlj] . While the coher- 
ence peak can be fully suppressed only for lar ge e lectron- 
phonon couplings in clean superconductors |13| , it has 
been shown that the peak can be further suppressed in 
disordered superconductors and is absent in the region of 
strong disorder for only moderate couplings |l4[ . 

The microscopic interplay of disorder, magnetic fluc- 
tuations and superconductivity is reflected in the behav- 
ior of the qp inelastic lifetime. In this paper we present 
a calculation for the qp inelastic scattering rate due to 
spin fluctuations within a formalism developed in previ- 
ous works The calculation is undertaken by first 
obtaining an effective fluctuation propagator in the su- 
perconducting state, and then using the exact eigenstate 
formalism as used in the case of Coulomb scattering with 
the replacement of the Coulomb propagator and vertex 
with the derived fluctuation propagators and vertices. It 
is shown that the rate is qualitatively similar to the rate 
due to Coulomb interactions with addition of the Stoner 
enhancement. Finally we discuss our results in terms of 
STM data on the borocarbide superconductors. 

II. CALCULATIONS 

The scattering rate from paramagnons in clean super- 
conductors on a lattice is well known for the case of s— 
or d— wave superconductors | p6| . The calculation for the 
inelastic scattering rate due to paramagnon exchange in 
disorder metals is also well known . In both cases the 
results are similar to the scattering rate from long-range 
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Coulomb interactions, with an additional enhancement 
via the Stoner factor 1/(1 — /), where / = UNp, U is 
a phenomenological short range interaction, and Np is 
the DOS per spin at the Fermi level. For dirty metals 
and superconductors, the electron-phonon interaction is 
reduced via "collision drag" relative to the clean case 
psf , while the electron-electron interaction is enhanced 
by disorder due to the breakdown of screening by diffusive 
electrons |l^,|l^ . The latter enhancement of the scatter- 
ing rate at the Fermi surface is p^/''^{Ep /T)^/'^ compared 
to that of 3D clean materials. Here p is the dimensionless 
measure of disorder, with p = pj pM-, where p is the ex- 
trapolated residual resistivity and pu the Mott number, 
which in a jellium model is given by pu — 'ii^'^ / e^kp. We 
use units such that fcs = 7i = 1. However, calculations 
for the scattering rate calculated for superconductors on 
a lattice have treated impurities and interactions in- 
dependently and therefore do not capture the disorder 
enhancement derived for conventional superconductors. 
Therefore in this paper we investigate the interplay of 
disorder, superconductivity, and magnetism by revisiting 
the problem of inelastic scattering. 

The spin fluctuation propagator is given by the sum of 
longitudinal K■\^^ and transverse K^^ paramagnons, re- 
spectively They can be expressed in terms of the 
polarization bubble x 



K^^ = UxU + UxUxK^^. 

Kn= u + UxKn- 



(1) 



Solving these equations we obtain the fluctuation propa- 
gator i(q, w) 



t{q,Lu) = if It + = 



i-c/V(q,^^) i-c/x(q,f^)' 

(2) 



However in the superconducting state one must distin- 
guish between charge and spin response couplings due to 
their different coherence factors. Therefore in the super- 
conducting state the propagator splits into two contribu- 
tions given by pl[ 



2 1-[/x.(q,f^) 



1 U^Xcjcj,^) 

2i + uxc{ci,n)' 



(3) 



with Xc,s the charge, spin susceptibilities, respectively. In 
the following we perform calculations in the continuum 
limit and neglect lattice effects. This is certainly impor- 
tant in order to capture strong scattering via qp exchange 
of antiferromagnetic reciprocal lattice vector momenta 
Q. However, the incipient magnetic instability via para- 
magnon exchange nevertheless is reflected via the Stoner 
criterion. Albeit a naive approach to the borocarbides 
or other materials with strong antiferromagnetic fluctu- 
ations, the results allow us to qualitatively estimate the 



effects of disorder on qp inelastic scattering from param- 
agnons. 

The gauge-invariant charge polarization Xc has been 
calculated in disordered superconductors in Ref. ||2^. It 
has the structure x(q, t^) = B{c[,uj) + Bc{ci,uj). Here B 
is the density response function in the pair approxima- 
tion , while Be contains the collective excitation (the 
Anderson-Bogolubov mode) which restores gauge invari- 
ance. It was shown that for fc^C S> 1 collective effects 
can be ignored and that the " pair approximation" for the 
polarization is adequate, where ^ = ^l/mnAp^^ is the 
dirty-limit coherence length. For T = the polarization 
can be written as 



X'M,^) = 0"(q, V^(c^ - 2A))e(c. - 2 A) 



{uj + 2A)E{a) 



UJ + 2A 



K{a) 



(4) 



while the spin susceptibility is given by the Mattis- 
Bardeen result IM 



X';(q,c.) = 0"(q, ^uj{u-2A))Q{u - 2 A) 

X [{uj + 2A.)E{a) - ^AK{a)] . (5) 

Here, a = , and E and K are complete elliptical 

integrals of the first and second kinds, respectively. 0" 
is the spectrum of the density Kubo function for non- 
interacting electrons. It can calculated by a variety of 
techniques for various limits of disorder. For clean met- 
als, the spectrum is white, 



(q, ej 



clean, 



(6) 



while for diffusive qp dynamics, 4>" is given by a diffusion 
pole 



'(q,e)-A^F 



diffusive, 



(7) 



with D the diffusion constant. Here we have neglected 
Cooper propagator renormalization, which can be shown 
to give a smaller contribution to the scattering rate 
than Diffusion propagator renormalization by a factor of 
l/kp^- _ ^ 

The limiting behavior for finite temperatures with 
T << A is given as: 



X^',,(q, n«A)^ f7</."(q, V2Mi)e-^/^ 

, charge. 



(A/T) ln(4r/r!), 



spin, 



Xc.s 



(q, n>2A)K A0"(q, 2V2A)^/^^TjAe-'^/'^ 
1/2, charge, 



1, 



spm. 



(8) 



(9) 



Thus the behavior of the spin and charge susceptibilities 
yields different contributions to the paramagnon scatter- 
ing in the charge and spin channel. 
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The paramagnon contribution to the self energy can 
be spht in the usual way into an anomalous and even 
and odd normal pieces. It has been shown for the case 
of long-range Coulomb interactions that the even part 
of the normal self energy contribution can be ignored, 
and can be shown for the spin-fluctuation case as well 
I psf . Expanding near the qp pole in the BCS Green's 
function [^5| , we obtain the expression for the on-shell 
inelastic scattering rate due to paramagnons exchange in 
the charge channel Tc and spin channel F^, 

q 

dx 

X / —[I{x + uj)+n{x)]t';^[ci,x) 
Jo I" 

X [G"(e, uj + x)± A/ujF"{e, uj + x)] + {uj ^ -cj), 

where n, f are Bose and Fermi distributions, respectively, 
G" and F" are the imaginary parts of the bare normal 
and anomalous BCS Green's functions, respectively, Z' 
is the real part of the qp renormalization, and —uj) 
denotes the addition of terms which differ from the ones 
written only by the sign of uj. 

Substituting Eqs.(3-9) into Eq.(lO), we obtain the in- 
elastic scattering rate = 2Fs. It can be shown that 
the contribution to the scattering rate from the charge 
channel yields a subdominant contribution for all val- 
ues of disorder and interaction / < 1 compared to the 
long-range Coulomb contribution calculated in Rcf . . 
Therefore for the remainder of the paper we neglect Fc 
and focus on Fs . The scattering rate is dominated by qp 
population at the gap edge. For T = 0, an injected qp 
must have enough energy to give up to break a Cooper 
pair (3 A) and for << 1 we obtain. 



rj=°{n > 3A)= 



3/^7r-^ A A A /n 
16(1- 7)3/2 Z'ii^B^ 



clean, 

(f)'''^/¥> dirty. 



(11) 



with F{x) = x{x'^/2 
2) ln[a; 



- l)v/(a;- 2)2-1 + {x/2 - 
1]. At finite temperatures, 



2 + -v/(^-2)2 
the Cooper pair recombination rate is dominated by the 
kinematic factor F-'^ cx e~2'^/^ and qp scattering rate 



cx e 



-A/T_ Pqj. a qp at the gap edge, the dominant 



contribution to the recombination rate is given by 



Ff (A >> T) 



37r2/- 



2 t2 



T2 



8a/2(1 - 1)3/2 Z'Ef 
7==, clean, 



-2A/T 



(12) 



(f ) dirty, 
while for the scattering rate we obtain to leading order 



Ff(A >> T) 



37r/2ln(2) A 
8(1-7)3/2 Z'V 2A' 
^ clean, 





-A/T 



dirty. 



(13) 



We see a similar behavior between the paramagnon and 
long-range Coulomb contributions to the inelastic scat- 
tering rate Jist . Tg possesses the same temperature de- 
pendence as the Coulomb contribution, with the expo- 
nential temperature dependence reflecting the necessity 
of two qps per scattering event. Further, we see the same 
disorder enhancement (p3/2(£'^/A)^/2 ) relative to the 
clean case as in the long-range Coulomb case. Lastly, we 
note that the energy gap A acts as a cut off for the di- 
vergence of the rate that occurs in the 2 — d dirty normal 
calculation just as in the long-range Coulomb case 
]ll. 

On top of the disorder enhancement, there is the 
Stoner enhancement relative to the Coulomb contribu- 
tions due to the nearness of a magnetic instability. In 
materials close to the instability, this contribution will 
be dominant over the Coulomb and phonon terms except 
for very low temperatures, where the power-law tem- 
perature dependence of the phonon contribution takes 
over ||l5|,^. We note that our expression are valid for 
A/Ep << 1 — 7 << 1, i.e., provided that one is not 
too close to the Stoner criterion for magnetism, 7=1. 
At the instability, the rate saturates as it does in the 
case of a normal metal near the metal-insulator transi- 
tion |2^. However in order to accurately describe the 
dynamics at the magnetic transition one needs to use a 
more sophisticated spin fluctuation propagator than the 
one derived here from RPA diagrams only, which tend to 
overestimate paramagnon effects p7[ |. 

Finally, we can compare the results to the values of 
the scattering rates inferred from STM data on clean 
and thin films of LuNi2B2C. To our knowledge, a tem- 
perature dependence of the scattering rate has not yet 
been published, nor has a reliable estimate of the scat- 
tering rate been made from optical (Raman or infrared) 
or Hall probes as has been done in the high Tc cuprates. 
Moreover no systematic study of the effects of impuri- 
ties and doping have been made concerning the scatter- 
ing rate. Nevertheless we can estimate if inelastic scat- 
tering from paramagnons in an s-wave superconductor 
is sufficient to explain the broadening observed in STM 
measurements Q . As a rough estimate for 1/ts we take 
7 ~ 2/3, Z' = 1/2, Fermi velocity vp ~ 3.5 x 10^ cm/s, 
Fermi energy Ep ^ 0.3 eV given from LDA estimates 
for LuNiaBzC from Ref. (2||. STM data taken at low 
temperatures in Ref. |Q| gives A = 18cm~^, which is 
consistent with Raman measurements jsj. This yields a 
scattering rate for clean systems at T = 0.5Tc from Eq. 
(13) of 1/ts = 1.3 X 10-3meV, or I/t^A = 6 x 10"^, 
which is clearly too small to match experiments. Either 
the scattering is most likely due to electron-phonon col- 



3 



lisions |g9[ or perhaps due to large gap anisotropy. 

Since p(T = 0) increases quickly as Co is doped in 
lH, rising by over an order of magnitude for 15% Co 
doping Ql , it may be feasible that the disorder enhance- 
ment for 1/ts in an s-wave scenario could lead to spec- 
tral weight at low frequencies observed via magnetic field 
anisotropy or Raman ||^ measurements. An estimate 
for the Mott number is difficult since the parameters 
vp^kp and the other parameters entering in Eqs. (12- 
13) are presumably disorder dependent, and it is not clear 
where the metal-insulator transition occurs for this com- 
pound. A conservative estimate from the loffe-Regel cri- 
terion in Ref. Q| gives pM ~ AOOfifl-cm. Therefore taking 
p{T = 0) - lOO/irj-cm as in Ref. [| into Eq. (13) only 
gives 1/tsA ^ 10^^, which is clearly too small to account 
for the large broadening observed via STM even in rel- 
atively clean films nor is it sufficient to account for the 
substantial spectral weight observed at low frequencies 
via Raman scattering. It is tempting to therefore con- 
clude that the large broadening comes either from nodal 
qps in conventional (extended s-) or unconventional (d-) 
pair states. 

However, there are problems in each scenario. Small 
amount of Co doping (on the few percent level) quickly 
push these materials into the dirty limit « 1) 0. 
If the gap possessed extended s— wave symmetry, the 
disorder would be sufficient to wash out any remaining 
anisotropy and necessarily lead to sharp threshold be- 
havior, which is not observed. On the other hand, the 
disorder would also lead to a sharp drop in Tc if the gap 
possessed d— wave symmetry and unconventional super- 
conductivity would be expected to be completely sup- 
pressed Q for 15% Co doping, which again is not ob- 
served. Therefore it is unclear from current data whether 
superconductivity is conventional or not, and perhaps the 
situation is clouded by the presence of additional non- 
superconducting bands, which would also yield a non- 
vanishing zero bias conductance and low frequency spec- 
tral weight. It would thus be extremely useful to study 
impurity and cation dopings further to determine if the 
enhanced scattering rates are responsible for the behav- 
ior indicative of unconventional pairing as the disorder 
is increased. Raman scattering measurements would be 
very useful in this regard, and remains a topic for further 
investigation. 
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